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We build a spherical halo model for galaxies using a general scalar-tensor theory of gravity in its 
Newtonian limit. The scalar field is described by a time-independent Klein-Gordon equation with 
a source that is coupled to the standard Poisson equation of Newtonian gravity. Our model, by 
construction, fits both the observed rotation velocities of stars in spirals and a typical luminosity 
profile. As a result, the form of the new Newtonian potential, the scalar field, and dark matter 
distribution in a galaxy are determined. Taking into account the constraints for the fundamental 
parameters of the theory (A, a), we analyze the infiuence of the scalar field in the dark matter 
distribution, resulting in shallow density profiles in galactic centers. 
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I. INTRODUCTION 



The recent indirect observational evidence of dark mat- 
ter (DM)and dark energy in the Universe [E H, H, 0, H, 
H, 0, H, Q has motivated the study of new cosmologi- 
cal and astrophysical scenarios that can encompass these 
observations. At a cosmological level, the quintessence 
scenario [l^, UH, provides fittings to the present ac- 
celerated expansion rate of the Universe and should be 
consistent with the other above-mentioned observations. 
Typically, quintessence introduces new fields that have 
their origin in theories that attempt to unify all forces of 
nature (strings, braneworlds). These unification schemes 
result in extensions to General Relativity (GR) which de- 
termine a new dynamics. Scalar-tensor theories (STT) 
of gravity are examples of effective theories that stem 
from such unifying schemes [l^, As one may sus- 

pect, in addition to cosmological consequences, such the- 
oretical extensions also predict local, astrophysical ef- 
fects. Traditionally, to understand the d:^namics of a 
galaxy, a DM profile has been introduced 15] or, alter- 
natively, a modification to the Newtonian gravitational 
law, e.g. via Yukawa couplings [l^, or even modifications 
to the Newtonian motion law, such as MOND [13, fisj . 
Within the first approach, several DM candidates have 
been proposed, including a scalar field (SF) as DM itself 
[igl [20| . The second approach has been used to obtain 
flat rotation curves in spirals via STT, without using DM 
halos. Adjusting rotation curve profiles using STT im- 
plies that the strength of the Yukawa coupling, a, has to 
be negative, leading to a phantomlike nonminimally cou- 
pled field. Moreover, the adjustment of different rotation 
curves of specific galaxies, points to different values for 
the range parameter of the theory, A pH , implying a mass 
spectrum for the fundamental theory at hand. The third 
approach solves flat rotation curves dynamics in spirals, 
but fails to fully understand cluster dynamics. This lat- 
ter approach will not be considered here. Our model is a 



combination of the first two approaches mentioned above: 
we use both a DM halo and a SF nonminimally coupled 
to GR, as long as these two elements could simultane- 
ously play a role in the galactic dynamics, and mitigate 
the constraints imposed on STT parameters, when using 
that theory alone. 

Following this line of thought, we have recently studied 
SiSlthe influence that STT have at galactic scales, see 
also |24| . We considered the Newtonian limit of STT and 
computed potential-density pairs of a spherical galaxy in 
which NFW's [25| and Dehnen's [261 density profiles were 
used. We also computed some other relevant observa- 
tional quantities (rotation curves, dispersions), by which 
we accounted for the influence of the STT scalar flelds in 
the Galaxy. Such influence is characterized by the two 
parameters (A, a) of STT, see below. In this work we use 
this formalism to build a galactic model that is, by con- 
struction, consistent with the measured rotation curves 
of stars and with some luminosity profile. As a result, 
the form of the Newtonian potential is exactly solved, 
and the SF and DM distribution in a galaxy are numeri- 
cally computed, and their specific features depend on the 
fundamental parameters (A, a) of the STT. For the val- 
ues of the parameter space analyzed, the resulting DM 
has shallow profiles near the center. 

Some other models such as the gravitational suppres- 
sion hypothesis |27[ have been put forward, in which a 
Yukawa term is added in the Newtonian potential. Re- 
cently, [2^ analyzed the rotation curves best flt in this 
model, in which a NFW profile is used. They con- 
cluded that this hypothesis does not fit several rota- 
tional curves of spirals and, hence, does not solve the 
core/cuspy problem of DM in the center of galaxies. Ar- 
guments in favor of a cusp like center can be found in 

Refs. m a in, m, m, however recently more ev- 
idence has emerged favoring a corelike galactic center 
[34, 35, 36, 37, 38, M, SO, El El, S, H SI • In H 
a way can be found to reconcile both approaches. 

The present work is organized as follows: in the next 
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section we give a brief description of the STT Newtonian 
approximation, where the SF background value is set to 
have the usual gravitational constant value at small dis- 
tances, r <C A. In Sec. IIII[ we build a galactic model 
by giving the rotation curves and baryon density profile, 
and in Sec. IIVI we analyze the influence of the pair of 
parameters (A, a) on the SF and DM density distribu- 
tions. Finally, in the last section we discuss the results 
and present our conclusions. 



II. THE NEWTONIAN APPROXIMATION OF 
STT 

A typical STT is given by the following Lagrangian 
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^;=(^))-2a2^;=<,)((0)l^;=<^)-2l-^.<^)) anda ^ 1/(3+ 
2a; (0)) \^={£)\u!{(t)) is a generalization of the Brans-Dicke 
parameter [47|. 

In the above expansion we have set the cosmological 
constant equal to zero since within galactic scales its in- 
fluence is negligible. This is because the average den- 
sity in a galaxy is much larger than a cosmological con- 
stant that is compatible with observations. Thus, we 
only consider the influence of luminous and dark matter. 
These matter components gravitate in accordance with 
the modified-Newtonian theory determined by Eqs. (j4|) 
and ([5|). The latter is a Klein-Gordon equation which 
contains the boson field of mass m, whose Compton wave- 
length (A = h/mc) implies a length scale for the new dy- 
namics. We shall assume this scale to be of the order of 
kiloparsecs. 

Note that Eq. ([4]) can be cast as a Poisson equation 



+LM{g,.) , (1) for ^(/ioo + ^) 



where g^^, is the metric, is a SF, and ijj{<t>) and V{4>) 
are arbitrary functions of it. C-Mig^iu) is the matter La- 
grangian. From Eq. ([T]) one obtains the gravity and SF 
equations. Thus, the gravitational equation is 
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The SF Klein-Gordon equation is 
(t)V' - 2V 1 



3 + 2w 



(3) 



where a prime (') denotes the derivative with respect to 
SF. 

In accordance with the Newtonian approximation, 
gravity and SF are weak. Then, we expect to have small 
deviations of the SF around the background field. As- 
suming also that the velocities of stars and DM particles 
are nonrelativistic, we perform the expansion of the field 
equations around the background quantities (</)) and ry^^ , 
i-e., g^iu = Vp.i' + h^_,u and 4 = {(j)) + S<p. 

The equations governing the weak energy (Newtonian) 
limit of STT are well known [22, 51, Isfl, [HI HI] and 
written here in physical units (ft, 7^ 1, c ^ 1) 
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i7rp--V'6<j,\ ,(4) 
(5) 



where the background value is chosen such as {(j)) — 
(1 + a)c^ /Gn, a choice that sets the effective Newto- 
nian constant to the one locally observed, see [i^, [5^ 
for a detailed discussion. Eqs. ([4]) and ([5|) represent 
the Newtonian limit of a set of STT that are distin- 
guished by the effective square mass = a:{{(f>)V^'^^^^ — 
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thus, the new Newtonian potential is now given by 



(6) 
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Solutions to these equations, the so-called potential- 
density pairs [54], were found for the NFW's and 
Dehnen's density profiles [HI and for axisymmetric sys- 
tems [2^. For point particles the solution is well known, 
see for instance [22, 50] , and with the choice of the above- 
mentioned background field, one has: 

Gn M 



$Ar = -- 



1 + a r 



(8) 



where M is the point particle mass producing the field. 
The strength of the new scalar force is given by a and 
its action range by A. For local scales, r <C A, deviations 
from the Newtonian theory are exponentially suppressed, 
and for r A the Newtonian constant diminishes (aug- 
ments) to Gjv/(1 + ct) for positive (negative) a. Re- 
cently, the effect of STT has been investigated in different 
cosmological scenarios in which variations of the Newto- 
nian constant are constrained. For instance, (ssj studied 
the infiuence of varying G^v on the Doppler peaks of the 
CMBR, and concluded that their parameter (^ = G/Gn) 
should be in the interval 0.75 < ^ < 1.74 in order to be 
within the error bars of the CMBR measurements. In our 
notation this translates into —0.43 < a < 0.33. However, 
this range for a has to be taken as a rough estimation, 
since these authors have only considered a variation of 
Gn, and not a full perturbation study within STT. The 
latter has been done by [56] , who found some deviations 
from the Newtonian dynamics, that when translated into 
our strength parameter would correspond to a = 0.04; 
however, they do not compare their results with obser- 
vations. On the other hand, a structure formation anal- 
ysis has been done by (STj . in which deviations of the 
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matter power spectrum are studied by adding a Yukawa 
potential to the Newtonian. They found some aUowed 
dynamics, that turn out to constrain our parameter to 
be within —1.0 < a < 0.5; but again a self-consistent 
perturbation study in general STT is missing. Thus, the 
above three estimates can be taken as order-of-magnitude 
constraints for our models. For definitiveness, we will 
take values within the range —0.3 < a < 3. The value 
a = —0.3 yields an asymptotic growing factor of 1.4 in 
Gtv, whereas the value a = 3.0 makes Gjv to asymptoti- 
cally reduce by one-fourth. 



implying that Mdm ~ 10~^^eV. The smallness of this 
mass would have to be explained by a particle physics 
theory, e.g. similar to light scalar presented by [60|, yet 
nonminimally coupled to gravity. 

Both baryons and DM 'feel' the same gravitational po- 
tential, <i>Ar, but are differently distributed in the galaxy. 
For the baryon component we assume a Freeman-disk 
density profile [UHl], that is, 



PB{r) = 



-r/rd 



(9) 



III. A GALACTIC MODEL 

We proceed to build a galactic model by assuming 
that the total matter content consists of two components, 
baryons and cold DM, Pt = Pb + Pdm! Baryons rep- 
resent stars, and the cold DM component could be of 
any type. The dynamics is determined by the theory 
explained in the preceding section. There are two pos- 
sibilities on the DM origin: i) DM is not related to the 
SF, and ii) DM is associated with the boson produced 
by the SF. In the former case, DM can be for example 
an ensemble of neutralinos, whose mass is in the range 
200GeV < < 300GeV [Hi], within an effective super- 
gravity theory that nonminimally couples to gravity, see 
for example [59|. In the latter case, ii), the mass of the 
DM particle is given by Mdm — h/\c, with A ~ kpc. 



where Md is the mass of the disk and its radius. 

For the DM density we do not assume a particular pro- 
file. Instead, we proceed to find its form by imposing a 
general rotation pattern. In the past attempts have been 
made to determine a universal rotation curve (URC) pro- 
file, beginning with the pioneering work of Ref. [63]. In 
the nineties the authors of Refs. J54, jj^l considered more 
than 1100 optical and radio data of Sb-Im spirals to find 
a phenomenological URC profile valid out the outermost 
radius where data were available at that time. Recently, 
they have considered more data and have modified the 
profile [6^, extending it out to its virial radius, that is, 
including the DM halo part. This profile is supposed to 
be valid for spirals of different types [67| but a number 
of issues are still open [6^ . Accordingly, we assume for 
our model that stars and DM particles obey the following 
URC profile [11]: 
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where the functions I and K are the modified Bessel func- 
tions, and po and ro are the scaling density and radius of 
the Burkert density profile [68j. The first part accounts 
for the disk contribution and the second for the halo's. In 
a previous work we have assumed a simpler, flat rotation 



curve profile [53|. 

The given circular velocity determines the form of 
gravitational potential ^n, through Eq. (jlOp . which in 
turn is related to pdni and Scfi through Eqs. (|4l5l7p . 

Integrating Eq. (fTU|) for $jv yields 
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which leads to the motion of test particles in the Galaxy. Substituting this result in the original system, Eqs. ^ [5]) 
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transform into the following two equations 

^'Sct>~-^^S4, = -^V2$^, (12) 
[l + a) Gn 

1 o 1 / rnc \ 2 ^ PnTn 1 / mc \ ^ ^ 



for two variables, p_dm and Scj). The second equality of 
Eq. (|13p results since the disk contribution of Laplacian 
of $Ar cancels out with baryon density, see Eq. (fTTj) . 
Thus, the resulting DM profile is the Burkert profile plus 
the SF contribution. By substituting Eq. pTjl into Eq. 
p2)) . one solves for the scalar perturbation. Then, us- 
ing Eq. (fT3|) one solves for the DM profile. The results 
depend on the STT parameters (A, a) and on the rota- 
tion curve fitting parameters {pQ,ro, Md,rd). Particular 
galaxies fix the latter parameters. As an example, we 
had chosen a set of four galaxies of different types that 
have been used to test particular gravity theorie s |28l| and 
to test the validity of some density profiles [H, l70j |. The 
galaxies are: DD047, ESO 116-G12, NGC 7339 and UGC 
4325. A greater galaxy set could have been used to test 
the rotation profile given by Eq. (jlOp. but this has al- 
ready been done in Ref. [Ig , and in Ref. [6^ additional 
arguments are given in favor of this rotation profile. In 
general, flat density profiles, such as Burkert's, tend to 
match rotation curves data within Newtonian dynamics 
[sl] . In table |T] we present the properties and best fitted 
values of the galaxies' parameters, and in Fig. [T]the best 
fits are shown. The rotation curve data fit quite well to 
the profile URC given by Eq. ^TU\i . In doing the fitting 
we have taken the values for rd, which is directly mea- 
sured by optical observations, given in Ref. [i^ for the 
first three galaxies and the values given in Ref. [70| for 
UGC 4325. Then, we have varied Md, ro and po- 

The first term of the right-hand-side of Eq. is 
the Burkert profile, pBUR = Po ''o/(^ + ro){r^ + r^) and 
one can define = g^^j^^^) (^) ^'l> that to express 
Eq. (Uni) as pDM = PBUR + P4>- The density p^ is the 
contribution of the SF fluctuation to the total density. 
Given this, the total density is px — pn + Pdm = Pb + 
Pbur ~^ Ptp- In what follows we proceed to numerically 
integrate the above equations. 



IV. SPHERICAL SOLUTION 

Unfortunately, these quantities cannot be computed 
analytically, since the Newtonian potential involves cylin- 
drical and spherical dependencies. Therefore to solve the 
Eq. p2p we will assume spherical symmetry. Then, we 
use the integrals given in Ref. [22] to numerically solve 
for 6(j). 

We perform the integration for the LSB galaxy UGC 
4325; because the NEW model does not properly fit with 



the observations, so we plan for this galaxy to contrast 
our results with them. 

Let us explain the range of values for the STT pa- 
rameters taken in our analysis. Originally, STT [i^, 13 
were thought for positives values of a to have a stan- 
dard kinetic term in Eq. ^j. But negatives values are 
also theoretically possible [Til [t^ and they have been 
applied, for instance, to accomplish, without a poten- 
tial in Eq. ([T]), an inflationary era in isotropic "TS] and 
anisotropic models [t^ . or more recently, to explain the 
present accelerated expansion of the Universe in some 
quintessence models [tsI. [76|. Thus, we will consider pos- 
itive as well negative values of a subject to the constraints 
mentioned at the end of section|lTl Therefore, we will an- 
alyze the solutions in the interval —0.3 < a < 3.0. On 
the other hand, we will assume that A is in the interval 
0.1 kpc < A < 50 kpc to fit galactic scales. 



A. Solutions for positive a 

In Fig. [2^ we plot the resulting density profiles for 
a = 3.0 and A = 1.0 kpc. Of particular interest is the 
form of the DM profile that flattens near the center of the 
Galaxy. The latter is again shown in Fig. [2}d, where for 
comparison the standard Newtonian density profile, pN, 
is plotted, which is the profile to have Vijj^c'H given by 
Eq. (fTU]) . but turning off the SF; i.e., pN is the density 
found by solving the Poisson equation with Newtonian 
potential given by ^nh in Eq. p4|) . Thus, the New- 
tonian density is the Burkert profile One observes 
that Pdm is always bigger than p^, since the effect of 
the SF is to diminish the effective gravitational constant 
for r > A, since GcS = G(l -I- ae~'~/^)/{l + a), thus to 
compensate the reduction of the gravitational constant, a 
denser DM profile is necessary to have the same rotation 
curve profile. Concerning the behavior of the profile, we 
have computed a numerical fit of the inner part of the 
curve (r <C rd), showing that it behaves approximately 
as PDM ~ r-'^""' with 7DM ~ 0.00006t;5:^g!^^i; the uncer- 
tainties stemming basically from the uncertainties in tq. 
On the other hand, the standard Newtonian model is just 
a httle steeper pN ~ r"^" with -fN « O.OOOlOlEJ JJgl^EJ^. 
Both of these profiles are essentially shallow. The NEW 
that best fitted rotation curves data is included for com- 
parison, which is known to be cuspier in the inner re- 
gion, pnfw ~ ^--1 00001 Pqj. j. ^ ^jjg behavior fol- 
lows PDM - r-^"" with 5dm ~ 0.040l°;°J^, pN ~ r-^" 
with 5n ~ Omitoml. The NEW model behaves as 
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TABLE I: Properties and best fitting parameters of the galaxies used. 



Galaxy 


Type 


rd [kpc] 


Ma [Me] 


ro [kpc] 


PO [Mg/kpC^] 


Xred 


DDO 17 


IB 


0..") 


:!.(){) = ().()2 X 10' 


.").l:! ± 0.09 


2.07 ± O.O:! X 10' 


1.74 


ES0116-G12 


SBcd 


1.7 


2.09 ± 0.08 X 10'-* 


4.77 ±0.03 


4.44 ± 0.04 X 10^ 


0.99 


NGC 7339 


SABd 


1.5 


1.10 ±0.01 X 10^° 


3.03 ± 0.03 


1.60 ± 0.02 X 10* 


1.69 


UGC4325 


SA 


1.6 


8.42 ± 0.47 X 10® 


40.54 ± 10.16 


6.59 ± 0.08 X 10^ 


3.56 




FIG. 1: The continuous line is the fitted rotation curve for data of the four galaxies. The short-dashed line is the exponential 
disk contribution and the long-dashed line is DM's. At the bottom we plotted the residuals {vohs — vurc)- 
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PNFW ^ r-''«^"' with 6nfw ^ 1-003. For r - tq the 
behavior follows pdm ^ r^^oM .^^jth 5dm ~ l-44to 26' 
- r-*" with (Jat « 1.44t°J^, and the NFW model 
behaves as pnfw ^ t^^nfw ^j^jj 5nfw ~ 1-08. Fig. 
[2}d shows that the DM profile is bigger than NFW's be- 
yond r ~ 0.6 kpc. On the other hand, Ref. [t^ shows 
that NFW fits are bad for this galaxy, because arbitrary 
low concentrations (i.e. large NFW scaling length r^) are 
needed. In our case, the best fitted curve implies a con- 
centration c — 1.22, which is clearly inconsistent with the 
cosmological expected values c > 5 (for a wide range of 
V200 values) [1, [t^- Thus, what we essentially see in our 
plots of the NFW profile is its cuspy region, as shown in 
Fig. Wp- 

The solution shown for the SF is the interior solution 
and eventually at some r = R the exterior solution is 
valid [131; one could think of R to be of the order of 
the halo size. Thus for r > R, the SF exponentially 
vanishes and we obtain a standard Newtonian behavior. 
Therefore, asymptotically, for r ^ rg, puu ^ similar 
to pnfw ~ T-^- 

In Figs. [2^-b we have plotted the SF and DM profiles 
for various A. The SF fluctuation diminishes going from 
the center to outer parts, and for smaller A the decay is 
stronger in inner regions. Fig. [3^ shows also that the SF 
is bigger for larger A, since for large A, GeS approaches to 
G. A consistency check implies that the SF must comply 
with the condition 6(j) < (I + a)c^Gj^ — (1 -|- a) 2.1 x 
lO^^Mo/kpc in order to validate the perturbation theory 
used, and it is indeed the case for the pair of parameters 
(A, a) chosen. On the other hand, the DM profiles in 
Fig. [Jb diminish by augmenting A. The DM profile more 
deviates from the standard Newtonian one for smaller A. 

In Figs. S^-b we plotted the SF and DM profiles for 
various a and fixed A = f .0 kpc. Again the constraint on 
5(t> is fulfilled (Fig. 2^). As expected, for small a the DM 
profile tends to the standard Newtonian one (Fig. \^). 

B. Solutions for negative a 

In Fig. [S^ we again plot all densities, as in Fig. [2^, but 
now for a — —O.I and A — I.O kpc. The DM profile in 
the inner regions (r ^ r^) is given by poM ~ 
with 7dm ~ 0-000Illg'QQoo2, that is quite similar as 
the standard Newtonian, which gives pN ^ r^''" , with 
7Ar ~ 0.000I0lo oo°o2; tliey are shallow. In contrast, 
NFW's profile in the same region that has a cuspy power 
^^pyi, = -1. 00001 . Fig. shows the DM, New- 

tonian and NFW profiles for comparison. It is clear 
that the DM profile is less massive than the standard 
Newtonian. This is because the effective gravitational 
function, Gcg — G(I -I- ae~^^'^)/{l + a), is increased 
for negative a, and the gravitational pull, being pro- 
portional to the term GeffM, is compensated by a de- 
crease in M(r), that is by a decrease of the DM profile. 
For r ^ Td the behavior is as follows: pdm ~ r^^oM 
with 5dm ~ O-OSSIq'oJq, the standard Newtonian pro- 



file is pn ~ r-^" with 5n « 0.0411° °^^ the NFW's 
behaves nearly as Snfw = 1-00. For r ^ ro the be- 
havior follows pdm ^ r^^DM .^^rith 5dm ~ l-44to;26' 
PN r-''" with 5n ~ l-44t[J;^^, the NFW model be- 
haves as pnfw r"''™^'*' with 5nfw ~ 1-08. As in 
the a-positive case, Fig. [SJd shows again that beyond 
some region {r ^ 1 kpc) the DM profile is bigger than 
NFW's. Asymptotically, for r ^ tq, beyond some point 
the exterior solution is valid, thus pdm ^ r^^, similar to 
Pnfw ^ r~^- 

Figs, in^-b show the SF and DM profiles for various 
A, respectively. The behaviors show systematic tenden- 
cies: the smaller A, the smaller \54>\, and again and for 
smaller A, the decay is stronger in inner regions (Fig. [5^). 
The SF complies with 5(j) < {1 + a)c^Gj^ , guaranteeing 
the validity of the perturbation approach. On the other 
hand, the bigger A, the closer the DM profile approaches 
to the standard Newtonian model (Fig. Wp). 

In Figs. [7^-b we plotted the SF and DM profiles for 
various negative a and fixed A = 1.0 kpc. Again the 
constraint on 54) is fulfilled. As expected, for smaller 
I a I the DM profile tends to the standard Newtonian one 
{a = 0). 



V. DISCUSSION AND CONCLUSIONS 

We have considered a general STT in its weak energy 
limit in which two free parameters appear (A, a). This 
pair can be constrained by the above-mentioned obser- 
vations: the first parameter is the Compton wavelength 
associated with a light boson particle, which we have 
taken to be of the order of kiloparsecs; this value implies 
a change of the Newtonian constant only at distances of 
the order of galactic scales or bigger, and therefore does 
not conflict with local deviations of Newtonian dynam- 
ics [zll. The other parameter is the strength of the new 
scalar force, given by a, which is subject to cosmologi- 
cal constraints. Accordingly, we have taken values for a 
within the range —0.3 < a < 3. We do not use bigger 
positive (negative) values for a because they predict a 
weaker (stronger) gravitational constant on scales larger 
than A, and bring some unacceptable cosmological effects 
[5i,[5i,[57|. 

Using the STT formalism, we have constructed a galac- 
tic model with a distribution of stars and DM that obey 
a rotation profile compatible with observations given by 
the URC of ^66^, as is shown in Fig. [T] This fitting 
determines the form of the effective gravitational poten- 
tial, $Ar, given by Eq. pT|) . Then, by taking a typical 
density profile for baryons (Freeman exponential profile) , 
the only two quantities to be determined are the SF fluc- 
tuation and DM profile, which are given by Eqs. (fT^ 
and (|f 3p . respectively, for which we have found numer- 
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FIG. 3: (a) It is shown the SF perturbation and (b) the DM profile for various A and fixed a = 3.0. 
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ical solutions'^. In a previous work [53l | we have found 
analytic solutions for a flat rotation curve profile. Those 
solutions are similar to NF W's [l^ . The numerical solu- 
tions found here are shallow near the galaxy center, since 
the URCH, Eq. stems from the Burkert profile [H], 

which is shallow for r ^ rp. In order to make a compar- 
ison with the standard Newtonian model, we turned off 
the SF and solved for the density in the standard Pois- 
son using Eq. (fTTj) . The resulting profile is called pN 
and is the well known Burkert profile. The DM profile 
found here is slightly shallower (cuspier) than the New- 
tonian profile for positive (negative) values of a, and in 
any case much shallower than NFW's. To quantitatively 
analyze the solutions we have considered the two scales 
involved in the URC fitting, that are the disc radius, 
Trf, and the DM scale, tq. The DM solutions near the 
galactic center (r <ti rd) with a positive are a bit shal- 
lower [p ^ r '^■'^'^"'J^-o.ooooi ) than those with a negative 

_p, ^nni 1 +0 00004 

[p ^ r '^■""'^^^-0.00002 )^ and both of them are much shal- 
lower than NFW's profiles [H, [H, HH [H . Solutions at 
r ^ rd for a positive decay as (pdm ^ r -o ooTj and 
for a negative decay weaker [p ^ r •'"-o ooo ) . At r ~ ro 
for positive a the DM profile behaves as p ~ r"' *'*-" 
and for negative a the behavior is p '--^ r ^■**-o.26. The 
uncertainties stemming from the allowed values of the 
fitting parameters reported in table I. For the allowed 
values of rg , the DM exponents for r <^ rd vary between 
20%-34%, for r^rd between 17%-29%, and for r > 
vary between 18% - 24%. The slopes of the DM profiles 
change smoothly, so that their behavior does not qualita- 
tively differ from their mean value. With respect to the 
variations of Md and po the exponents vary less than (or 
of the order of) 1%. 

On the other hand, for positive (negative) a, poM is 
bigger (smaller) than p^ always, since the effect of the 
SF is to diminish (augment) the effective gravitational 

^ In other models, such as the gravitational suppression hypothesis 
|27( one assumes a NFW DM profile to determine the parameters 
(A, a) of that theory by fitting theoretical curves to the rotation 
curves data [Ssll . In our approach we determine the DM profile 



constant for r > A, being GcS = G{1 + ae~^^^)/{l + a), 
thus it is necessary to compensate it with a corresponding 
larger (smaller) DM density. Finally, asymptotically, for 
r ^ ro, the exterior solution is valid, and thus pdm ^ 
is similar to pnfw ~ 

We have found numerical solutions for the allowed pa- 
rameter space for strength of the SF potential, —0.3 < 
a < 3.0, and within galactic distances, 0.1 kpc < A < 50 
kpc. The results indicate some systematic tendencies: 
the effect of the SF is more apparent for A < tq, and its 
influence attenuates for A > rg, since in this regime the 
behavior is essentially Newtonian (A = oo) for r < tq. 
For small |a| the DM profile tends to the standard New- 
tonian one (a = 0). 

The intention of the present work was to study the 
influence of a massive SF in a galaxy model that is com- 
patible with a typical baryon distribution and follows the 
URC of observed galaxies. This construction fixes the 
Newtonian potential, and for the STT parameters ana- 
lyzed, the resulting DM profile is shallow at the center 
of the halo. These results are encouraging showing the 
important role that both contributions, DM and STT, 
could have in the dynamics of the systems. We think that 
numerical N-body simulations using the STT of gravity 
have to be done to confirm the solutions discussed here. 
Some preliminary computations have been carried out in 
Ref. [7g|,l80|. 
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